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Some rigorous results can be derived using a very simple approach to hadron spec- 
troscopy, in which a static potential is associated with non-relativistic kinematics. Several 
regularities of the experimental spectrum are explained by such models. It is underlined 
that certain methods developed for hadronic physics have found applications in other 
fields, in particular atomic physics. A few results can be extended to cases involving 
spin-dependent forces or relativistic kinematics. 

1. Introduction 

The dynamics of light quarks is presumably very intricate, with non-perturbative effects 
even at short distances, and highly-relativistic motion of the constituents inside hadrons. 
Nevertheless, it is interesting to consider a fictitious world, where the quark dynamics 
is governed by a simple Hamiltonian with non-relativistic kinematics supplemented by a 
static, flavour-independent potential. Most regularities derived from the properties of the 
Schrodinger equation are observed in the experimental spectrum. This suggests that the 
QCD theory of quark confinement should lead to similar regularities. 

Several results derived in the context of quark models have been successfully applied to 
other few-body problems, for instance in atomic physics. Examples will be given in this 
review. 

A challenge consists of extending the theorems on level order, convexity, etc., to less 
naive Hamiltonians with spin-dependent forces and relativistic kinematics. Some of the 
first results will be mentioned. 

2. Results on mesons 

Studies of the quark model have been stimulated by the discovery of Jf$ and T res- 
onances and their excitations and, in particular, the simultaneous description of both 
spectra using the same potential (flavour independence). Many rigorous results have 
been summarized in the reviews by Quigg and Rosner |I| and by Martin and Grosse [^,[3] . 
A few examples are given in this section, dealing with energy levels. A result on the wave 
function will be mentioned in Sec. |9|. 
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Current potentials reproduce the observed pattern that E(1P) < E(2S). Within the 
notation (n, £) adopted here for quarkonium, the radial wave function has n — 1 nodes 
and the principal quantum number of atomic physics is n + £. It has been proved that 
E(n + l,£) > E(n, i + 1) if AV > 0, and the reverse if AV < 0. The Coulomb degeneracy 
is recovered as a limiting case. From the Gauss theorem, the sign of the Laplacian 
AV = r~ 2 [r 2 V'}' reflects whether the charge Q(r) seen at distance r grows (asymptotic 
freedom), decreases or remains constant. 

This "Coulomb theorem" can be applied successfully to muonic atoms, which are sen- 
sitive to the size of the nucleus {Q{r) /*), and to alkaline atoms whose last electron 
penetrates the inner electron shells (Q(r) \). The latter property also hold for the 
metastable p — He molecules described by Korobov at this conference M. 

Another theorem explains how the harmonic oscillator (h.o.) degeneracy E(n + 1,£) = 
E(n, £ + 2) is broken. A strict inequality is obtained if the sign of V" is constant ||. 

In both the complete Hamiltonian p 2 / (2m 1 ) + p 2 ,/ (2m 2 ) + V^(r 12 ) or its reduced version 
p 2 /(2/x) + V(r), the individual inverse masses, raf 1 , or the inverse reduced mass, fi^ 1 , 
enter linearly through a positive operator p 2 . This implies that each energy level is an 
increasing function of this inverse mass, and that the ground-state energy (or the sum of 
first levels) is a concave function of this variable || . Let us mention two among the many 
applications. 

The dynamics of heavy quarkonia QQ is very sensitive to the change from Q = c to 
Q = b, while open-flavour mesons are dominated by the light quark mass. Thus B — D is 
essentially the difference between quark masses h — c, while the effect of this difference is 
reduced in T — \I/ by the change in the binding energy. After removing or estimating the 
hyperfine effects, one gets from the data 

-(T - tf) ~ 3.2 GeV < (B- D) ~ 3.3 GeV. (1) 

The second example is concerned with comparing symmetric mesons QQ with hidden 
flavour to asymmetric mesons Q'Q with open flavours. As (be), for instance, has an 
inverse reduced mass which is the average between that of (bb) and that of (cc), then 
[0)0 (be) > (bb + cc)/2 is expected. Moreover, d-E/d/i -1 = (p 2 )/2, and this latter matrix 
element can be estimated from the excitation spectrum and bounded along its path from 
(bb) or (cc) to (be). One eventually ends with a rigorous window for the mass of the 
spin-averaged ground state |3||| 

6.26 < (6c) < 6.43 GeV. (2) 

Of course, hyperfine effects should be added in more precise tests. The accuracy of the 
first measurement 

(6c) = 6.4 ± 0.4 GeV (3) 

is not yet sufficient to test flavour independence. 

Convexity in the inverse reduced mass also implies that (bq 2 ) + (cqi) < (<?i) + (cg 2 ) if 
b > c and g 2 > Qi- In particular, 



(6c) > (6s) + (cc) - (cs) ~ 6.52 GeV, 



(4) 
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which is less efficient but much more easily obtained than the r.h.s. of Eq. fl2|). 

Note that an inequality such as D 2 + H 2 < 2HD should hold for variants of the hydrogen 
molecule involving isotopes, as they experience the same Born-Oppenheimer potential. 

One of the "sides" of the Coulomb theorem can be adapted to the Klein-Gordon equa- 
tion @. If AF < 0, then E(n+1,£) < E(n,£ + 1). 

3. Level order of baryon spectra 

For many years, the only widespread knowledge of the 3-body problem was the harmonic 
oscillator. This remains true outside the few-body community. The discussion on baryon 
excitations is thus often restricted to situations where V = v { r ij)i with v(r) = Kr 2 + 5v, 
and Sv treated as a correction. 

First-order perturbation theory is usually excellent, especially if the oscillator strength 
K is variationally adjusted to minimise the magnitude of the corrections. However, when 
first-order perturbation is shown (or claimed) to produce a crossing of levels, one is rea- 
sonably worried about higher-order terms, and a more rigorous treatment of the energy 
spectrum becomes desirable. 

In this context, it should be recalled that in the N and A sector (not in the A one!), 
the radial excitation with positive parity (Roper resonance) comes first, and the orbital 
excitation with positive parity just above. On the other hand, any simple potential such 
as v (r) = Xr — a/r gives the opposite ordering. This problem has been known for many 
years. It has been vigorously revisited by Glozman and collaborators |nj and many 
others. The true nature of light-quark dynamics is outside the scope of this review, which 
will be focused on the level order in 3-body potential models. 

A decomposition better than V — Yl K r h + $ v is provided by the generalised partial- 
wave expansion 

V = V (p)+5V, (5) 

where p oc (r 2 2 + r| 3 + r 2 ^) 1 / 2 is the hyperradius. The last term SV gives a very small 
correction to the first levels |ITJ . With the hyperscalar potential Vo only the wave function 
reads ^ = p~ 5//2 w(p)P[£](f2), where the last factor contains the "grand-angular" part. The 
energy and the hyperradial part are governed by 

U "(P) ~ + m[E - Vo(p)]u(p) = 0, (6) 

very similar to the usual radial equations of the 2-body problem, except that the effective 
angular momentum is now i = 3/2 for the ground-state and its radial excitations and 
£ = 5/2 for the first orbital excitation with negative parity. The Coulomb theorem holds 
for non-integer I. If AV > 0, then E(2S) > E(1P), i.e., the Roper comes above the 



orbital excitation |12j| . Note that a three-body potential cannot be distinguished from a 
simple pairwise interaction once it is reduced to its hyperscalar component Vq by suitable 
angular integration. It also results from numerical tests that relativistic kinematics does 
not change significantly the relative magnitude of orbital vs. radial excitation energies. 

The problem of the Roper resonance is thus a difficult one if one requires a model to 
reproduce perfectly the observed masses. Several solutions have been proposed, including 
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accepting with philosophy that a 100 MeV discrepancy on a particular level does not 
necessarily means the end of the conventional quark level ||13|| . One might mention: 
1. Introducing in the potential some terms with negative Laplacian. In particular, scalar- 
meson exchange induces a Yukawa-like potential v oc — exp(— fir)/r. Its strength is how- 
ever not sufficient to change the level order, when it is associated with a linear potential 
of reasonable magnitude 



2. Introducing a specific spin or spin-flavour dependence in the potential, with the con- 
sequence that Eq. (||) is no longer valid: the effective potential in negative-parity states 
differs from the one governing the ground state and its radial excitations. This is what 
occurs with pion exchange in the models introduced by Glozman and others. 

3. More drastic solutions consist of considering the Roper resonance as a multiquark state 
or a hybrid state, possibly mixed with conventional radial excitations. The production 
and decay properties do not provide much support for such an interpretation. 

Back to spectral theory, the splitting of levels in the nearly hyperscalar potential (|J) 
is very similar to the famous pattern of the N = 2 h.o multiplet flTRITIJ, except that 



the Roper is disentangled fl^ , |TT]j . A similar result is found for higher negative-parity 
excitation: the split N = 3 levels of the nearly harmonic model are separated into a 
radially excited L = 1 and a set of split L = 3 levels. Details and references can be found 
in Refs. HJT^jn . 



4. Test of flavour independence for baryons 

A look at the baryon spectrum gives guide lines to our intuition which are contradictory. 
We have just seen that the "vertical spectrum", i.e., the ordering of radial vs. orbital 
excitations suggests flavour-dependent potentials such as those mediated by exchange of 
Goldstone bosons. On the other hand, the "horizontal spectrum", i.e., the spectrum of 
ground states with different flavour content, is rather smooth. There is a very continuous 
and slow increase in mass when a strange quark replaces an ordinary one. Strangeness 
excitation usually costs less than spin excitation. 

This latter observation suggests describing all ground states with the same universal 
(flavour independent) potential and then applying moderate spin corrections "a la DGG" 
[pSfl, where flavour dependence occurs only through the mass of the constituent in the 



Breit-Fermi term. This was basically the substance of the Isgur-Karl model and 
other attempts with non-harmonic confinement p0 |. 



This leads to considering the mathematical properties of the 3-body spectrum in a given 
potential as a function of the constituent masses. We have seen in Sec. |2] that for any 
2-body potential V which is flavour independent (just required to produce bound states), 
inequalities such as (bb) + (cc) < 2 (be) can be derived. The analogue for baryons, 

(QQq) + (Q'Q'q) < 2(QQ'q), (7) 

requires, however, mild restrictions, which are always satisfied by the potentials used in 
quark model, whatever quark masses are used for Q, Q' and q. Flavour independence is 
the only serious hypothesis. See Ref. [|1J . 

For instance, the equal spacing rule, Or — E* = S* — S* = S* — A, is understood as 
follows: the central force gives a concave behaviour, with for instance Q~ — 2* < H* — £*, 
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but a quasi perfect linearity is restored by the spin-spin interaction which acts more 
strongly on light quarks. A similar scenario holds for the Gell-Mann-Okubo formula [20 . 

Inequalities can also be written for baryons with heavy flavour, some of them being 
more accessible than others to experimental checks in the near future. Examples are 

3{bcs) > (bbb) + (ccc) + (sss), (8) 
2(bcq) > {bbq) + (ccq), 2(cgg) > (ccq) + (qqq)- 



5. Baryons with two heavy flavours 

There is a renewed interest in this subject [p2p3l| . The recent observation of the (be) 
mesons demonstrates our ability to reconstruct hadrons with two heavy quarks from their 
decay products. 

Baryons with two heavy quarks (QQ'q) are rather fascinating: they combine the adi- 
abatic motion of two heavy quarks as in J/\l/ and T mesons with the highly relativistic 
motion of a light quark as in flavoured mesons D or B. 

The wave function of (QQq) exhibits a clear diquark clustering with r(QQ) <C r(Qq) 
for the average distances. This does not necessarily mean that for a given potential model, 
a naive two-step calculation is justified. Here I mean: estimate first the (QQ) mass using 
the direct potential v(QQ) only, and then solve the \{QQ) — q] 2-body problem using a 
point-like diquark. If v is harmonic, one would miss a factor 3/2 in the effective spring 
constant of the (QQ) system, and thus a factor (3/2) 1//2 in its excitation energy. 

On the other hand, it has been checked that the Born-Oppenheimer approximation 
works extremely well for these (QQq) systems, even when the quark mass ratio Q/q is 
not very large. This system is the analogue of in atomic physics. 



6. S^r-expansion and stability of multiquarks 

The deuteron and other nuclei are bound by long-range nuclear forces, not by direct 
quark-quark interaction, and thus cannot be considered as genuine multiquark states. 
Stable multiquarks have been searched for in all possible channels, with the exception 
of very exotic hadronic configurations with several heavy flavours, and the results are 
negative. 

The absence, or at least the non-proliferation of stable multiquarks, is reproduced in 
most model calculations, even the simplest ones. 



I shall not discuss too much here the multiquark states such as the hexaquark H [[M] or 
the pentaquark P [[25[] whose binding is tentatively achieved by chromomagnetic forces. 
For these states, the predictions vary from one author to another, and the questions are: 
i) Is the one-gluon-exchange appropriate for describing hyperfine effects ? Note that the 
Goldstone-boson-exchange model gives slightly different predictions on multiquarks p6 



as compared to the chromomagnetic model, but the main message remains identical: there 
is no proliferation of multiquarks in the spectrum. 

ii) What is the value of the short-range correlation in a multiquark, as compared to 
ordinary hadrons? This quantity governs, indeed, the strength of chromomagnetic effects. 
Another issue deals with the absence of collective binding due to chromoelectric (con- 
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fining) forces. Consider a four-body problem with Hamiltonian 

4 2 

Ha{x) = J2 7T- + C 1 - 2a; )(^ + + (1 + x)(V 13 + V u + V 23 + V u ) 
i=i Zm 

= H S + xH MS , (9) 

where the parameter x measures the departure from a fully symmetric interaction. From 
the variational principle, the ground-state energy E(x) is maximal at x = 0. Indeed, with 
obvious notation, 

E(x) = (V(x)\H(x)\V(x)) < (tf(0)|#(a;)|tf(0)) = E(0), (10) 

where the last step results from \P(0) being symmetric and thus insensitive to the mixed- 
symmetry part, Hms- As x enters H linearly, E(x) is a concave function and thus 
should exhibit near x = a parabola-like shape. This means that 

E(0) > E(xx) > E(x 2 ) if 0<X!<a;2 (11) 

is guaranteed, and 

E(x[) > E(x 2 ) if x[ < < x 2 and \x[\ < x 2 (12) 

is very likely. In short, a larger asymmetry leads to a lower algebraic energy. 

Consider now the most primitive ansatz for the multiquark interaction, the so-called 
additive colour model 

v = -^E^M^) = >>,'•('•<,)• (is) 

i<j i<j 

This is a pairwise interaction with a colour dependence that corresponds to pure octet 
exchange. The normalisation is such that v is the central potential for quarkonia. If 
this model is applied to four-quark configurations (qqqq) assuming a frozen colour wave 
function, then when one looks at the distribution of the strength factors gij, one observes 
that the largest asymmetry is obtained for the meson-meson threshold, for which an 
appropriate labelling gives 

= 0, except 5-12 = g 3i = I. (14) 

It is easily checked that the strength distributions corresponding to a "true baryonium" , 
with (qq) in a 3 colour state, and (qq) in 3 colour state, or a "mock baryonium" state, 
with a 6-6 colour structure, are more clustered around the symmetry point where all 
are equal to 1/3. (Colour neutrality ensures that J29ij is the same for all colour-singlet 
configurations.) 

Indeed, solving the equal-mass (qqqq) problem using a frozen colour wave function and 
the potential ( [13] ) never leads to stable multiquarks. 
Binding requires a combination of 

— mixing of several colour wave functions 

— spin-dependent effects 

— a better estimate of the multiquark potential. Note however that the recent lattice 
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estimate by Green et al. |27| suggests that the true 4-quark potential is less attractive 
than the ansatz ([13]). 

— another type of Sn breaking which would be favourable to multiquark binding. 

This latter mechanism was suggested several years ago. Its most remarkable feature is 
that theorists agree on the main result, namely that {qqQQ) becomes stable in the limit 
where the quark mass ratio Q/q becomes very large. 

The Hamiltonian 

ff = ^(pi+ri) + ^(ri+ri) + v, (15) 

with V even under charge conjugation, i.e., simultaneous (1 <-> 3) and (2 «-> 4) exchanges, 
can be seen as the x = 1 value of 

H(x) = H cvcn + xH odd , 



\4m 4M 



2 2 2 

P 2 - P 3 - P4 



When x grows from to 1, the meson-meson threshold is unchanged, because the inverse 
reduced mass of (m, M) is already an average of m _1 and M _1 . On the other hand, the 
four-quark energy E{x) decreases. Stability of {qqQQ) becomes possible if Q 3> q. 
Remarks on this multiquark are in order: 

1. If the mass ratio approaches the critical value for stability, the wave function has large 
meson-meson components. For larger Q/q, the binding is more pronounced, and the 
system acquires another structure in space and colour. The two heavy quarks form a 
colour 3 source (as in every Sqq baryon), which in turn, form a colour singlet with the 
two antiquarks, as in every flavoured antibaryon (Qqq). One thus uses well established 
colour coupling. This contrasts with the highly speculative colour structure in "colour- 
chemistry" . 

2. The same mechanism holds for (+,+,—,—) configurations in atomic physics. The 
positronium molecule Ps2 is very weakly bound. The hydrogen molecule, where different 
masses experience the same potential as Ps2, is much more stable. In fact the stability of 
Ps2, as first demonstrated in 1947 by Hylleraas and Ore implies stability for hydrogen 



and any similar configurations such as K + K + n~n~ |29|. On the other hand, the molecule 
ppe + e~ is unstable, even in the limit where strong interaction and electromagnetic anni- 
hilation is switched off: one can perform an Sn expansion analogous to (|IBD; the 4-body 
energy decreases as asymmetry (m _1 ^ M _1 ) is implemented, but the threshold benefits 
much more from this effect. In the limit where M 3> m, the system breaks in a point-like 
protonium which cannot attract the remaining positronium. The critical mass ratio for 
stability has been estimated to about M/m ~ 2.2 |BU|. 



7. Hall— Post inequalities for mesons and baryons 

The simple additive model ( |T3"D gives for baryons the "1/2" rule 

Vqqq = 12 V q<l( r ij), with V qq = ^V q g, (17) 

i<3 
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or, for Hamiltonians, after distributing the kinetic energy among different terms 



Hqqq ^ ^11 2 ^ 

i<j i<j 



P 2 P 2 

— + — + V-(r--' 
2m + 2m + qqKl3J 



The variational principle, using the baryon wave function gives 

E(qqq) > 3E(qq)/2, (19) 
or for the masses, including the contribution of the constituents, 

(qqq)/3 > (qq)/2. (20) 

In other words, matter is heavier per quark in baryons than in mesons. This is seen, 
e.g., by comparing fi~(1672) and 0(1020), where we are dealing with the spin-triplet 
potential, including the spin-spin term. Otherwise, one can compare the spin 1/2 baryon 
to a suitable combination of spin and spin 1 mesons. 

Some generalisations to unequal quark masses are straightforward, for instance 

2(bcs) > (be) + (cs) + (sb) (21) 

On the other hand, if one rewrites (|20|) as 

{qqq) + (qqq) > 3(g?), (22) 

which expresses the fact that quark rearrangement is energetically allowed for annihilation 
at rest, one observes that for a large enough mass ratio Q/q, the inequality becomes 
inverted, namely 



(QQQ) + (qqq) < 3(Qq). (23) 

For instance, a triple-charm antibaryon would not "annihilate" on ordinary matter. 

These inequalities were rediscovered in quark model studies [j3lH 33H but they were known 
from previous works on the stability of matter and the thermodynamic limit of iV-body 
systems p4|, and studies of few-nucleon systems |j35 |. 

The simplest result is the following. If 

N 2 

H N (m,g) = Y,^ + 9Y,n^) (24) 

i zm i<j 

governs a system of N bosons and En denotes the ground-state energy, then the identity 
H N (m, g) = ^ £ Ff_, (^m, g) , (25) 
leads to 

N /N — 1 \ 

E N (m, g) > Jj~zr^ E N-i ( -y^m, 9 J ■ (26) 
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This inequality is never saturated, because the overall translation energy of the (N — 1)- 
body subsystems within the iV-body system is neglected. An improvement consists of 

by the identity 

N 



replacing the decomposition 

1 



\m 



H N (m,g) 



(i) 
N- 



N-l 



m,g 



N-2 

relating the translation-invariant Hamiltonians 

Ep*] 2 



H N = H 



N 



2Nm 

This leads to the new inequality 

N 



E N (m,g) > 



N-2' 



-E 



N-l 



N 



N-l 



m,g 



(27) 



(2* 



(29) 



As Nm/(N - 1) is smaller than (N - l)m/(N - 2) for iV > 3, the bound (gp| is better 
than (p6|) , since any binding energy in a given potential is a decreasing function of the 
constituent mass. 

By recursion, an inequality such as (2B) or (|2l 



provides a bound on E^ in terms of the 
2-body energy Ei. In the case of the harmonic oscillator for bosons, the bound derived 
from fl2"9|) is saturated, whereas the bound derived from ( |26| ) is, by a factor \/2, smaller 
than the exact result at large N. 

The generalisation to different constituent masses has been carried out [^j. The de- 
composition of the Hamiltonian involves some free parameters which can be adjusted to 
optimise the inequality. Again, saturation is obtained in the case of harmonic forces. 



The situation is far more difficult for fermions. In Ref. p4 |, the decomposition 

„2 



H N (m,g) 



2?£ 



(N — l)m 



+ gV(r jt 



(30) 



is written down, which expresses Hn in terms of Hamiltonians with (N — 1) independent 
particles. For fermions, one obtains 



N 

E N (m,g) > —f N . 



{m(N-l),g), 



(31) 



where /jy is the cumulated energy of a system of N independent fermions, a notation 
borrowed from Ref. f38|. In Eq. (|3l|) , the translation energy of the (iV — l)-body subsystem 
is neglected. 

Basdevant and Martin [38| succeeded in improving this inequality in the case of power- 
law potentials V oc r n with n > 1 or combinations of those with positive weights. They 
obtained saturation for any number N of fermions in the harmonic limit n = 2. Recently 
Juillet et al. |$9[ used group-theoretical considerations to analyse the structure of iV-body 
wave function in terms of the permutation of a (N — l)-subset and derived the following 
inequality 

«?/ x N- 1 

E N(m, g) > Ar/Ar oWoo , 7T x (32) 



N(N -2)(2S +1 
S(N + 2S + 2)E S N _\ /2 



m, 



Ng 



N 

where S is the spin of the iV-fermion system. 



T ) + (S+1)(N-2S)E 



5+1/2 
N-l 



m . 



Ng 



N 
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8. Convergence of variational calculations 

When computing the binding energy of a baryon or a more complicated system in a 
given potential model, one is faced with the choice of a method for solving the few-body 
problem. 

Almost all possible methods have been applied to baryons. We already mentioned the 
approximation of a perturbed oscillator. The hyperspherical expansion and the Faddeev 
equations have also been applied very convincingly. The approximation of Feshbach and 
Rubinow [^(J , which looks at the best function of the perimeter r ij , is n °t very accurate 



for baryons with linear confinement, while it turns out to be very powerful for weakly- 
bound states with short-range interaction. 

Many variational methods are base on a decomposition 

M 

* M = ]Tc^(a, /?,...), (33) 

i=l 

where the linear parameters q are determined from a eigenvalue equation and the range 
parameters a,/3, . . . adjusted by suitable minimisation of the energy. For most choices of 
the basis functions </?j, one achieves a decent convergence of the variational energy E 
toward the exact energy E, with the consequence that the approximate variational wave 
function overlaps better and better the exact one, i.e., []5| 

as E M — > E. (34) 

However, the variational wave function \1/ M might differ appreciably from the exact \I/ in 
regions of negligible contribution to the normalisation integral but of interest for physics. 
In particular, the value of the wave function for — > governs the production rate, the 
electromagnetic and annihilation width of quarkonia, the contribution of VT-exchange to 
the weak decay of flavoured mesons and baryons, the parity-violating effects in atoms, 
etc. A comparison of the estimates of the correlation coefficient 

5 12 = {*\8®(r 12 )V) (35) 

using different 3-body methods is rather interesting. It shows that: 

1. The quality of the Faddeev wave function is impressive, even with a small number of 
partial waves included. 

2. The hyperspherical formalism provides a slow but safe convergence for this quantity. 

3. Some variational expansions of type (^) fail in reproducing 5 12 ! This is observed 
in particular when the basis is built out of the successive eigenstates of the harmonic 
oscillator. (These wave functions fall off too rapidly at large separation. To keep the 
normalisation and the average energy correct, a distortion is introduced at small separa- 
tion). 4. However, the expansion on a basis of Gaussians does not suffer from the same 
difficulties, as contributions of longer and longer range can always be introduced in the 
expansion. 

9. Bounds on short-range correlations in baryons 

As seen in the last section, it is rather difficult to calculate the 3-body (or higher) wave 
function with good zero-range properties. 
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To avoid lengthy computations, one can estimate 812 using a trick proposed by Schwinger 
for the 2-body problem |I] and generalised to larger systems [fil~| , fl2"|| . The basic idea is 
that A and Ag = A + [H, B] have same expectation value within any eigenstate of H . A 
suitable choice of B might make {^! m \Ab\^ m ) converge faster than (fy M \A\^ M ) to the 
exact (\&|yl| 1 I r ) as the size M of the basis is increased. 

For a 2-body problem with a central potential V(r), the operator A = 8^(r) is replaced 
by Ab oc dV/dr, which is much less sensitive to a local flaw of the wave function. In 
particular 8± 2 comes out exact for a linear potential, however bad is the trial wave function, 
provided it is normalised. One also understands why in such a purely linear potential, 8™ 2 
is the same for all ^-states, independent of the radial number n. A better result has been 
established 0: 5™ 2 increases with n if V" > 0, and decreases if V" < 0. For instance, 
current potentials —a/r + br explain why the leptonic coupling is smaller for radially 
excited quarkonia than for the ground state. 

For baryons, As contains also a small centrifugal term, besides the derivative of the 
full potential V with respect to r\2- Cohen and Lipkin |4"3"] have compared the £ and H 
hyperfme splittings in models where this is induced by a contact term. In other words, 
they compared 812 in (qqs) and (ssq). On can also 03] derive bounds on <5 12 . 

The short-range correlation between electrons in ordinary atoms, the annihilation prob- 
ability of exotic molecules with positrons can also be calculated using the Schwinger rule. 
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